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EXPRESSION MONITORING FOR HUMAN 
CYTOMEGALOVIRUS (HCMV) INFECTION 

BACKGROUND OF THE INVENTION 

5 Many biological functions are accomplished by altering the expression of 

various genes through transcriptional (e.g. through control of initiation, provision 
of RNA precursors, RNA processing, eta) and/or translational control. For 
example, fundamental biological processes such as cell cycle, cell differentiation 
and cell death, are often characterized by the variations in the expression levels of 

10 groups of genes. 

Gene expression is also associated with pathogenesis. For example, the 
lack of sufficient expression of functional tumor suppressor genes and/or the over 
expression of oncogene/protooncogenes could lead to tumorigenesis (Marshall, 
Ce//,64: 313-326(1991); Weinberg, Science, 254: 1138-1 146 (1991), incorporated 

15 herein by reference for all purposes). Thus, changes in the expression levels of 

particular genes (e.g. oncogenes or tumor suppressors) serve as signposts for the 
presence and progression of various diseases. 

The study of gene expression in the art has been generally concentrated on 
the regulatory regions of the gene of interest and on the relationships among a few 

20 genes. A number of transcriptional factors/DNA binding proteins have been 

identified and a limited number of regulatory pathways have been discovered. 
However, the expression of a particular gene is frequently regulated by the 
expression of a large number of other genes. The expression of those regulatory 
genes may also be under the control of additional genes. This complex regulatory 

25 relationship among genes constitutes a genetic network. The function and 

regulation of a particular gene can be best understood in the context of this genetic 
network. As the Human Genome Project and commercial genome research 
progress at a great rate, most, if not all, of the expressed genes will be partially 
sequenced in the near future. Understanding the functions and regulatory 

30 relationships among the large number of genes is becoming a difficult task with 
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traditional tools. Therefore, there is a need in the art to develop a systematic 
approach to understand the complex regulatory relationships among large numbers 
of genes. 



Summary of the Invention 



5 This invention provides methods, compositions, and apparatus for studying 

the complex regulatory relationships among host genes and viruses, in particular 
HCMV. In some of its specific applications, this invention provides methods, 
compositions, and apparatus for identifying drugs for preventing or ameliorating 
disease symptoms caused by HCMV. In other applications the invention provides 
10 methods for determining the stage of infection or the extent of tissue damage 

caused by HCMV infection. In another embodiment the invention provides a 
general method for narrowing large sets of genes which may be important down 
to smaller subsets of genes which have elevated probabilities of being biologically, 
physiologically, and medically relevant. 

15 Brief Description of the Drawings and Table 



Figure 1. Characterization of RNA target samples and reproducibility of 
array-based hybridization results. (Figure hi) 82, 68 and 51 probe pairs were used 
to interrogate the 5', middle and 3* portions of the GAPDH mRNA, which is 
constitutively expressed in fibroblasts. (Figure IB, Figure 1 Q The plots compare 

20 the average difference intensities (Avg. Diff. Intensities) of the 20 probe pairs 

interrogating each of the genes present in two independent experiments performed 
on the mock-infected cells (Figure 1A) or cells at 8 h after infection (Figure 1. E). 
The parallel lines flanking the center diagonal line indicate 3, 10 and 30-fold 
changes in intensity. With the exception of the thombospondin-1 gene in the 

25 mock-infected control, all other genes demonstrated an average difference in their 

hybridization intensities of less than 3-fold. 
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Figure 2A, Figure 2B, Figure 2C. Global survey of the differences in mRNA 
levels after HCMV infection. 

The plots show the variation in expression levels (Avg. Diff Intensities) between 
mock-infected cells and cells at 40 min, 8 h and 24 h after infection ( Figure 2A, 
5 Figure 2B, Figure 2C, respectively). Changes in expression of 3, 10, and 30-fold 

are highlighted by the parallel lines flanking the center diagonal line. 
Figure 3. Representative Northern blot analyses confirming changes in mRNA 
levels predicted by DNA array assay. Cultures of primary human diploid 
fibroblasts were infected with HCMV strain AD 169 or Toledo, and total cellular 
10 RNA was analyzed by Northern blot at 40 min, 8 h and 24 h after infection. Genes 

to which the probes correspond are identified to the right of the autoradiograms. 
M, mock-infected cells. 

Figure 4A and 4B (Table 1), Cellular mRNAs whose levels change by a factor of 
four or more after infection with HCMV. Identity of columns from left to right: 
1 5 GenBank accession number, name of gene encoding mRNA; time(s) after infection 

when a change in mRNA level was observed plus fold change; increase (U) or 
decrease (D) in steady state level of RNA; gene chip results confirmed in this 
report by northern blot (1), confirmed by another literature report (2), not 
confirmed (3). 
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DETAILED DESCRIPTION OF THE INVENTION 
I. Definitions 

Bind(s) substantially: "Bind(s) substantially" refers to complementary 
5 hybridization between a probe nucleic acid and a target nucleic acid and embraces 

minor mismatches that can be accommodated by reducing the stringency of the 
hybridization media to achieve the desired detection of the target polynucleotide 
sequence. 

Background: The terms "background" or "background signal intensity" refer to 

10 hybridization signals resulting from non-specific binding, or other interactions, 

between the labeled target nucleic acids and components of the oligonucleotide 
array (e.g., the oligonucleotide probes, control probes, the array substrate, etc.). 
Background signals may also be produced by intrinsic fluorescence of the array 
components themselves. A single background signal can be calculated for the 

15 entire array, or a different background signal may be calculated for each target 

nucleic acid. In a preferred embodiment, background is calculated as the average 
hybridization signal intensity for the lowest 5% to 10% of the probes in the array, 
or, where a different background signal is calculated for each target gene, for the 
lowest 5% to 10% of the probes for each gene. Of course, one of skill in the art 

20 will appreciate that where the probes to a particular gene hybridize well and thus 

appear to be specifically binding to a target sequence, they should not be used in 
a background signal calculation. Alternatively, background may be calculated as 
the average hybridization signal intensity produced by hybridization to probes that 
are not complementary to any sequence found in the sample {e.g. probes directed 

25 to nucleic acids of the opposite sense or to genes not found in the sample such as 

bacterial genes where the sample is mammalian nucleic acids). Background can 
also be calculated as the average signal intensity produced by regions of the array 
that lack any probes at all. 
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Cis-acting : The term "cis-acting" is used here to refer to the regulation of gene 
expression by a DNA subsequence in the same DNA molecule as the target gene. 
Cis-acting can be exerted either by the binding of trans-acting transcriptional 
factors or by long range control. 

5 Complexity: The term "complexity" is used here according to standard meaning 

of this term as established by Britten et al Methods of Enzymol. 29:363 (1974). 
See, also Cantor and Schimmel Biophysical Chemistry: Part III at 1228-1230 for 
further explanation of nucleic acid complexity. 

Hybridizing specifically to: The phrase "hybridizing specifically to" refers to 
10 the binding, duplexing, or hybridizing of a molecule substantially to or only to a 

particular nucleotide sequence or sequences under stringent conditions when that 
sequence is present in a complex mixture (e.g., total cellular) DNA or RNA. 

Introns: noncoding DNA sequences which separate neighboring coding regions. 
During gene transcription, introns, like exons, are transcribed into RNA but are 
1 5 subsequently removed by RNA splicing. 

Massive Parallel Screening: The phrase "massively parallel screening" refers to 
the simultaneous screening of at least about 100, preferably about 1000, more 
preferably about 10,000 and most preferably about 1,000,000 different nucleic acid 
hybridizations. 

20 Mismatch control: The term "mismatch control" or "mismatch probe" refer to a 

probe whose sequence is deliberately selected not to be perfectly complementary 
to a particular target sequence. For each mismatch (MM) control in a high-density 
array there typically exists a corresponding perfect match (PM) probe that is 
perfectly complementary to the same particular target sequence. The mismatch 

25 may comprise one or more bases. While the mismatch(s) may be located anywhere 

in the mismatch probe, terminal mismatches are less desirable as a terminal 
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mismatch is less likely to prevent hybridization of the target sequence. In a 
particularly preferred embodiment, the mismatch is located at or near the center of 
the probe such that the mismatch is most likely to destabilize the duplex with the 
target sequence under the test hybridization conditions. 

5 rpRNA or transcript: The term "mRNA" refers to transcripts of a gene. Transcripts 

are RNA including, for example, mature messenger RNA ready for translation, 
products of various stages of transcript processing. Transcript processing may 
include splicing, editing and degradation. 



Nucleic Acid: The terms "nucleic acid" or "nucleic acid molecule" refer to a 
1 0 deoxyribonucleotide or ribonucleotide polymer in either single-or double-stranded 

form, and unless otherwise limited, would encompass analogs of natural nucleotide 
that can function in a similar manner as naturally occurring nucleotide. A n 
oligo-nucleotide is a single-stranded nucleic acid of 2 to n bases, where n may be 
greater than 500 to 1000. Nucleic acids may be cloned or synthesized using any 
15 technique known in the art. They may also include non-naturally occurring 

nucleotide analogs, such as those which are modified to improve hybridization and 
peptide nucleic acids. 

Nucleic acid encoding a regulatory molecule: The regulatory molecule may be 
DNA, RNA or protein. Thus for example DNA sites which bind protein or other 
20 nucleic acid molecules are included within the class of regulatory molecules 

encoded by a nucleic acid. 

Perfect match probe: The term "perfect match probe" refers to a probe that has 
a sequence that is perfectly complementary to a particular target sequence. The test 
probe is typically perfectly complementary to a portion (subsequence) of the target 
25 sequence. The perfect match (PM) probe can be a "test probe", a "normalization 

control" probe, an expression level control probe and the like. A perfect match 
control or perfect match probe is, however, distinguished from a "mismatch 
control" or "mismatch probe." 
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Probe: As used herein a "probe" is defined as a nucleic acid, capable of binding 
to a target nucleic acid of complementary sequence through one or more types of 
chemical bonds, usually through complementary base pairing, usually through 
hydrogen bond formation As used herein, a probe may include natural (z.e. A, G, 
5 U, C, or T) or modified bases (7-deazaguanosine, inosine, eta). In addition, the 

bases in probes may be joined by a linkage other than a phosphodiester bond, so 
long as it does not interfere with hybridization. Thus, probes may be peptide 
nucleic acids in which the constituent bases are joined by peptide bonds rather than 
phosphodiester linkages. 

10 

Target nucleic acid: The term "target nucleic acid" refers to a nucleic acid (often 
derived from a biological sample), to which the probe is designed to specifically 
hybridize. It is either the presence or absence of the target nucleic acid that is to 
be detected, or the amount of the target nucleic acid that is to be quantified. The 
15 target nucleic acid has a sequence that is complementary to the nucleic acid 

sequence of the corresponding probe directed to the target. The term target nucleic 
acid may refer to the specific subsequence of a larger nucleic acid to which the 
probe is directed or to the overall sequence (e.g., gene or mRNA) whose expression 
level it is desired to detect The difference in usage will be apparent from context 



20 Trans-acting: The term "trans-acting" refers to regulation of gene expression by 

a product that is encoded by a gene at a remote location, usually as a result of 
binding to a cis-element. 



Stringent conditions: The term "stringent conditions" refers to conditions under 
which a probe will hybridize to its target subsequence, but with only insubstantial 
25 hybridization to other sequences or to other sequences such that the difference may 

be identified. Stringent conditions are sequence-dependent and will be different 
in different circumstances. Longer sequences hybridize specifically at higher 
temperatures. Generally, stringent conditions are selected to be about 5 0 C lower 
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than the thermal melting point (Tm) for the specific sequence at a defined ionic 
strength and pH. 

Subsequence: "Subsequence" refers to a sequence of nucleic acids that comprise 
a part of a longer sequence of nucleic acids. 

5 Thermal melting point (Tm): The Tm is the temperature, under defined ionic 

strength, pH, and nucleic acid concentration, at which 50% of the probes 
complementary to the target sequence hybridize to the target sequence at 
equilibrium. As the target sequences are generally present in excess, at Tm, 50% 
of the probes are occupied at equilibrium). Typically, stringent conditions will be 
10 those in which the salt concentration is at least about 0.01 to 1.0 M Na ion 

concentration (or other salts) at pH 7.0 to 8.3 and the temperature is at least about 
30 C for short probes (e.g, 10 to 50 nucleotide). Stringent conditions may also be 
achieved with the addition of destabilizing agents such as formamide. 

Quantifying : The term "quantifying" when used in the context of quantifying 
15 transcription levels of a gene can refer to absolute or to relative quantification. 

Absolute quantification may be accomplished by inclusion of known 
concentration^) of one or more target nucleic acids (e.g. control nucleic acids such 
as Bio B or with known amounts the target nucleic acids themselves) and 
referencing the hybridization intensity of unknowns with the known target nucleic 
20 acids (e.g. through generation of a standard curve). Alternatively, relative 

quantification can be accomplished by comparison of hybridization signals 
between two or more genes, or between two or more treatments to quantify the 
changes in hybridization intensity and, by implication, transcription level. 

25 Sequence identity: The percentage of sequence identity" or "sequence identity" 

is determined by comparing two optimally aligned sequences or subsequences 
over a comparison window or span, wherein the portion of the polynucleotide 
sequence in the comparison window may optionally comprise additions or 
deletions (i.e., gaps) as compared to the reference sequence (which does not 
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comprise additions or deletions) for optimal alignment of the two sequences. The 
percentage is calculated by determining the number of positions at which the 
identical subunit {e.g. nucleic acid base or amino acid residue) occurs in both 
sequences to yield the number of matched positions, dividing the number of 
5 matched positions by the total number of positions in the window of comparison 

and multiplying the result by 100 to yield the percentage of sequence identity. 
Percentage sequence identity when calculated using the programs GAP or 
BESTFTT (see below) is calculated using default gap weights. 
Methods of alignment of sequences for comparison are well known in the art. 

10 Optimal alignment of sequences for comparison may be conducted by the local 

homology algorithm of Smith and Waterman, Adv. Appl Math. 2: 482 (1981), by 
the homology alignment algorithm of Needleman and Wunsch J. Mol Biol 48: 
443 (1970), by the search for similarity method of Pearson and Lipman, Proa Natl 
Acad. Set. USA 85: 2444 (1988), by computerized implementations of these 

1 5 algorithms (including, but not limited to CLUSTAL in the PC/Gene program by 

Intelligenetics, Moutain View, California, GAP, BESTFIT, FASTA, and TFASTA 
in the Wisconsin Genetics Software Package, Genetics Computer Group (GCG), 
575 Science Dr., Madison, Wisconsin, USA), or by inspection. In particular, 
methods for aligning sequences using the CLUSTAL program are well described 

20 by Higgins and Sharp in Gene, 73: 237-244 (1988) and in CABIOS 5: 151-153 

(1989)). 

Up-stream or down-stream gene. If the expression of a first gene is regulated by 
a second gene, the second gene is called an "up-stream gene" for the first gene 
and the first gene is the "down-stream" gene of the second gene. The regulation 

25 of the first gene by second gene could be through trans-activation. For example, 

the first gene encodes a transcriptional factor that controls the expression of the 
second gene. The regulation can also be exerted by cis-acting. For example, 
the first gene is in the proximity of the second gene and exerts a positional 
effect on the expression of the second gene. In this case, the first gene does not 

30 have to be expressed in order to have an influence on the second gene. 
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According to the present invention, the stage of disease caused by HCMV 
infection can be determined. Expression levels of one or more genes which are 
induced or repressed by HCMV are determined in a first human cell sample. 
5 The first human cell sample comprises cells of a patient infected with HCMV 

and consists essentially of HCMV-infected cells. The expression levels of the 
one or more genes correlates with stage of disease progression of the HCMV 
infection. 

According to another embodiment of the invention the extent of tissue damage 

1 0 caused by HCMV infection is determined. Expression levels of one or more 

genes which are induced or repressed by HCMV are determined in a first human 
cell sample. The first human cell sample comprises cells of a patient infected 
with HCMV and consists essentially of HCMV-infected cells. The expression 
levels of the one or more genes correlates with extent of tissue damage caused 

15 by the HCMV infection. 

According to yet another aspect of the invention candidate drugs for 
preventing or ameliorating disease symptoms caused by HCMV are identified. 
Human cells are contacted with HCMV and a test agent. The contacting of the 
test agent and the HCMV with the cells can be at the same time or sequentially. 

20 Expression levels of one or more genes which are induced or repressed by 

HCMV are determined. Test agents are identified as candidate drugs if the test 
agents cause the human cells to express the one or more genes at the same(z.£, 
within 1 0-50% ) of the level at which the human cells express the one or more 
genes in the absence of HCMV infection. 

25 The genes whose expression levels are tested are those which are induced or 

repressed by HCMV. These are preferably those which are induced or repressed 
to a level which is at least two-fold, four-fold, eight-fold, ten-fold, fifteen-fold 
different than the level of expression in the absence of HCMV. More preferably 
the genes are selected from those genes identified in Table 1 as repressed or 

30 induced by HCMV. Of those genes which are identified in Table 1 , HLA-E, 

Ro/SSA, lipocortin-1, cPLA2, COX-2, thrombospondin-1 , mdMFTF, are 
preferred. 
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According to another aspect of the invention a general paradigm for 
identifying subsets of biologically, physiologically, or medically relevant genes 
is provided. This paradigm permits the prioritization of attention, investigation, 
and research on those genes which are most likely to have biologically, 
5 physiologically, or medically relevance. The paradigm involves the 

combination of expression data with other types of information contained within 
databases, including the general scientific literature, the patent literature, 
nucleotide and amino acid sequence databases, etc. By systematically searching 
databases for information which may shed light on the source of the expression 

10 data, one can make new connections and correlations which will confirm or 

suggest heightened biological, physiological, or medical relevance. The 
discussion section of the Example provides multiple examples of the operation 
of this paradigm. Generally, expression levels of at least two genes are 
compared between two cell samples. The two cell samples are the same but for 

15 a selected environmental, genetic, or developmental difference. These include 

without limitation contact with a drug or other exogenous chemical agent; 
temperature difference; mutation; viral infection; developmental stage 
difference; and bacterial infection. A set of genes whose expression levels 
differ between the two cell samples is identified. A database is searched to 

20 identify an environmental agent, gene, disease, biological phenomenon, or 

developmental stage previously associated with expression or loss of expression 
of individual members of the set of genes. As exemplified below, these may 
include without limitation an immunological reaction, a biochemical pathway, 
knock-out experimental animals, mutant animals or cells, diseases. When a 

25 common biological feature is identified between the selected environmental, 

genetic, or developmental difference used between the two sets of cells and the 
unselected environmental agent, gene, disease, or developmental stage 
identified from the database, the common gene is identified as being a member 
of a subset of genes which are improved targets for drug development. A 

30 common biological feature may be, for example, an association with a symptom 

of a disease, a phenotype, a reaction, etc, A computer-readable medium having 
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computer-executable instructions may also be used for performing the steps of 
this paradigm. Such a medium may be stored on disk or other suitable medium, 
as is well known in the art. 
The method by which expression levels are determined is not critical to the 
5 invention. Either mRNA or protein expression from one or more genes may be 

determined. Any method known the art for determining such expression levels 
can be used. These include without limitation hybridization to an array of 
oligonucleotides, serial analysis of gene expression, hybridization on a solid 
support, and immunological assays, such as Western blot, ELISA, 

10 immunohistochemistry. mRNA can be collected from the human cells, reverse 

transcribed, and used as a template for amplification. Any detection means 
known in the art can be used, including but not limited to use of a radioactive 
label, a fluorescent label, a chromophoric label, an enzymatic label. 
Various techniques known in the art render the screening of large numbers of 

15 genes'expression relatively straight-forward. Thus although a single gene's 

expression can be determined and can provide diagnostic and prognostic 
information, multiple genes can also be tested for their expression levels. In 
some embodiments at least 2,5,10, 15, 25, 50, 100, 150, 200, 350, 500, or 
1000 genes are tested to determine their expression levels. 

20 Cell samples for use in the assays of the present invention can be any cell type 

which is infected by HCMV. These include but are not limited to human cells 
fibroblasts, lymphocytes, epithelial cells, lung epithelial cells, and neuronal 
cells. When comparisons are done between an HCMV-infected cell sample and 
an uninfected cell sample, preferably the two samples are of the cell type. 

25 However this is not always possible. Comparisons can also be done with 

consensus expression profiles determined using a population of cells from 
normal, uninfected cells. 

Once expression levels of one or more genes are determined in a clinical 
human cell sample from HCMV-infected cells, they can be compared to the 

30 level found in uninfected controls, whether from the same individual and 

preferably the same type of cell, or from other individuals. A level of induction 
or repression greater than the predetermined threshold can be correlated with a 

12 



WO 00/011218 



PCT/US99/18772 



stage of disease progression or extent of tissue damage. Such correlation can be 
made by reference to a standard curve generated of stage of disease progression 
or extent of tissue damage as a function of gene expression levels. 
The methods of the invention involve quantifying the level of expression of a 
5 large number of genes. In some preferred embodiments, a high density 

oligonucleotide array is used to hybridize with a target nucleic acid sample to 
detect the expression level of a large number of genes, preferably more than 10, 
more preferably more than 100, and most preferably more than 1000 genes. 
Activity of a gene is reflected by the activity of its product(s): the proteins or 

10 other molecules encoded by the gene. Those product molecules perform 

biological functions. Directly measuring the activity of a gene product is, 
however, often difficult for certain genes. Instead, the immunological activities 
or the amount of the final product(s) or its peptide processing intermediates are 
determined as a measurement of the gene activity. More frequently, the 

1 5 amount or activity of intermediates, such as transcripts, RNA processing 

intermediates, or mature mRNAs are detected as a measurement of gene 
activity. 

In many cases, the form and function of the final produces) of a gene is 
unknown. In those cases, the activity of a gene is measured conveniently by the 
20 amount or activity of transcript(s), RNA processing intermediate^), mature 

mRNA(s) or its protein product(s) or functional activity of its protein 
product(s). 

Any methods that measure the activity of a gene are useful for at least some 
embodiments of this invention. For example, traditional Northern blotting and 
25 hybridization, nuclease protection, RT- PCR and differential display have been 

used for detecting gene activity. Those methods are useful for some 
embodiments of the invention. However, this invention is most useful in 
conjunction with methods for detecting the expression of a large number of 
genes. 

30 High density arrays are particularly useful for monitoring the expression 

control at the transcriptional, RNA processing and degradation level. The 
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fabrication and application of high density arrays in gene expression monitoring 
have been disclosed previously in, for example, WO 97/10365, WO 92/10588, 
U.S. Application Ser. No. 08/772,376 filed December 23, 1996; serial number 
08/529,1 15 filed on September 15, 1995; serial number 08/168,904 filed 
5 December 15, 1993; serial number 07/624,114 filed on December 6, 1990, serial 

number 07/362,901 filed June 7, 1990, all incorporated herein for all purposed 
by reference. La some embodiment using high density arrays, high density 
oligonucleotide arrays are synthesized using methods such as the Very Large 
Scale Immobilized Polymer Synthesis (VLSDPS) disclosed in U.S. Pat. No. 

10 5 ,445,934 incorporated herein for all purposes by reference. Each 

oligonucleotide occupies a known location on a substrate. A nucleic acid target 
sample is hybridized with a high density array of oligonucleotides and then the 
amount of target nucleic acids hybridized to each probe in the array is 
quantified. One preferred quantifying method is to use confocal 

15 microscope and fluorescent labels. The GeneChip® system (Affymetrix, Santa 

Clara, CA) is particularly suitable for quantifying the hybridization; however, it 
will be apparent to those of skill in the art that any similar systems or other 
effectively equivalent detection methods can also be used. 
High density arrays are suitable for quantifying a small variations in 

20 expression levels of a gene in the presence of a large population of 

heterogeneous nucleic acids. Such high density arrays can be fabricated either 
by de novo synthesis on a substrate or by spotting or transporting nucleic acid 
sequences onto specific locations of substrate. Nucleic acids are purified and/or 
isolated from biological materials, such as a bacterial plasmid containing a 

25 cloned segment of sequence of interest. Suitable nucleic acids are also 

produced by amplification of templates. As a nonlimiting illustration, 
polymerase chain reaction, and/or in vitro transcription, are suitable nucleic 
acid amplification methods. 
Synthesized oligonucleotide arrays are particularly preferred for this invention. 

30 Oligonucleotide arrays have numerous advantages, as opposed to other 

methods, such as efficiency of production, reduced intra- and inter array 
variability, increased information content and high signal-to-noise ratio. 
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Preferred high density arrays for gene function identification and genetic 
network mapping comprise greater than about 100, preferably greater than about 
1000, more preferably greater than about 16,000 and most preferably greater 
than 65,000 or 250,000 or even greater than about 1,000,000 different 
5 oligonucleotide probes, preferably in less than 1 cm 2 of surface area. The 

oligonucleotide probes range from about 5 to about 50 or about 500 
nucleotides, more preferably from about 10 to about 40 nucleotide and most 
preferably from about 15 to about 40 nucleotides in length. 

Massive Parallel Gene Expression Monitoring 

10 One preferred method for massive parallel gene expression monitoring is 

based upon high density nucleic acid arrays. Nucleic acid array methods for 
monitoring gene expression are disclosed and discussed in detail in PCT 
Application WO 092.10588 (published on June 25, 1992), all incorporated 
herein by reference for all purposes. 

1 5 Generally those methods of monitoring gene expression involve (a) 

providing a pool of target nucleic acids comprising RNA transcript(s) of one or 
more target gene(s), or nucleic acids derived from the RNA transcript(s); (b) 
hybridizing the nucleic acid sample to a high density array of probes and (c) 
detecting the hybridized nucleic acids and calculating a relative and/or absolute 

20 expression (transcription, RNA processing or degradation) level. 

(A). Providing a Nucleic Acid Sample 
One of skill in the art will appreciate that it is desirable to have nucleic 
samples containing target nucleic acid sequences that reflect the transcripts of 
interest. Therefore, suitable nucleic acid samples may contain transcripts of 
25 interest. Suitable nucleic acid samples, however, may contain nucleic acids 

derived from the transcripts of interest. As used herein, a nucleic acid derived 
from a transcript refers to a nucleic acid for whose synthesis the mRNA 
transcript or a subsequence thereof has ultimately served as a template. Thus, a 
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cDNA reverse transcribed from a transcript, an RNA transcribed from that 
cDNA, a DNA amplified from the cDNA, an RNA transcribed from the 
amplified DNA, etc., are all derived from the transcript and detection of such 
derived products is indicative of the presence and/or abundance of the original 
transcript in a sample. Thus, suitable samples include, but are not limited to, 
transcripts of the gene or genes, cDNA reverse transcribed from the transcript, 
cRNA transcribed from the cDNA, DNA amplified from the genes, RNA 
transcribed from amplified DNA, and the like. 

Transcripts, as used herein, may include, but not limited to pre-mRNA 
nascent transcript(s), transcript processing intermediates, mature mRNA(s) and 
degradation products. It is not necessary to monitor all types of transcripts to 
practice this invention. For example, one may choose to practice the invention 
to measure the mature mRNA levels only. 

In one embodiment, such sample is a homogenate of cells or tissues or 
other biological samples. Preferably, such sample is a total RNA preparation of 
a biological sample. More preferably in some embodiments, such a nucleic acid 
sample is the total mRNA isolated from a biological sample. Those of skill in the 
art will appreciate that the total mRNA prepared with most methods includes not 
only the mature mRNA, but also the RNA processing intermediates and nascent 
pre-mRNA transcripts. For example, total mRNA purified with a poly (dT) 
column contains RNA molecules with poly (A) tails. Those polyA + RNA 
molecules could be mature mRNA, RNA processing intermediates, nascent 
transcripts or degradation intermediates. 

Biological samples may be of any biological tissue or fluid or cells from 
any organism. Frequently the sample will be a "clinical sample" which is a sample 
derived from a patient. Clinical samples provide a rich source of information 
regarding the various states of genetic network or gene expression. Some 
embodiments of the invention are employed to detect mutations and to identify the 
phenotype of mutations. Such embodiments have extensive applications in clinical 
diagnostics and clinical studies. Typical clinical samples include, but are not 
limited to, sputum, blood, blood cells (e.g., white cells), tissue or fine needle 
biopsy samples, urine, peritoneal fluid, and pleural fluid, or cells therefrom. 
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Biological samples may also include sections of tissues, such as frozen sections or 
formalin fixed sections taken for histological purposes. 

Another typical source of biological samples are cell cultures where gene 
expression states can be manipulated to explore the relationship among genes. In 
5 one aspect of the invention, methods are provided to generate biological samples 

reflecting a wide variety of states of the genetic network. 

One of skill in the art would appreciate that it is desirable to inhibit or 
destroy RNase present in homogenates before homogenates can be used for 
hybridization. Methods of inhibiting or destroying nucleases are well known in the 

10 art. In some preferred embodiments, cells or tissues are homogenized in the 

presence of chaotropic agents to inhibit nuclease. In some other embodiments, 
RNase is inhibited or destroyed by heat treatment followed by proteinase treatment. 

Methods of isolating total mRNA are also well known to those of skill in 
the art. For example, methods of isolation and purification of nucleic acids are 

15 described in detail in Chapter 3 of Laboratory Techniques in Biochemistry and 

Molecular Biology: Hybridization With Nucleic Acid Probes, Parti. Theory and 
Nucleic Acid Preparation, P. Tijssen, ed. Elsevier, N.Y. (1993) and Chapter 3 of 
Laboratory Techniques in Biochemistry and Molecular Biology: Hybridization 
With Nucleic Acid Probes, Part L Theory and Nucleic Acid Preparation, P. 

20 Tijssen, ed. Elsevier, N.Y. (1993)). 

Li a preferred embodiment, the total RNA is isolated from a given sample 
using, for example, an acid guamdinium-phenol-chloroform extraction method and 
polyA 4 " mRNA is isolated by oligo(dT) column chromatography or by using (dT) 
on magnetic beads (see, e.g, Sambrook et al 9 Molecular Cloning: A Laboratory 

25 Manual (2nd ed.), Vols. 1-3, Cold Spring Harbor Laboratory, (1989), or Current 

Protocols in Molecular Biology, F. Ausubel et aL, ed. Greene Publishing and 
Wiley-Interscience, New York (1987)) . 

Frequently, it is desirable to amplify the nucleic acid sample prior to 
hybridization. One of skill in the art will appreciate that whatever amplification 

30 method is used, if a quantitative result is desired, care must be taken to use a 
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method that maintains or controls for the relative frequencies of the amplified 
nucleic acids to achieve quantitative amplification. 

Methods of "quantitative" amplification are well known to those of skill 
in the art. For example, quantitative PCR involves simultaneously co-amplifying 
a known quantity of a control sequence using the same primers. This provides an 
internal standard that may be used to calibrate the PCR reaction. The high density 
array may then include probes specific to the internal standard for quantification 
of the amplified nucleic acid. 

One preferred internal standard is a synthetic AW106 cRNA. The AW106 
cRNA is combined with RNA isolated from the sample according to standard 
techniques known to those of skilled in the art. The RNA is then reverse 
transcribed using a reverse transcriptase to provide copy DNA. The cDNA 
sequences are then amplified (e.g. 9 by PCR) using labeled primers. The 
amplification products are separated, typically by electrophoresis, and the amount 
of radioactivity (proportional to the amount of amplified product) is determined. 
The amount of mRNA in the sample is then calculated by comparison with the 
signal produced by the known AW106 RNA standard. Detailed protocols for 
quantitative PCR are provided in PCR Protocols, A Guide to Methods and 
Applications, Innis et aL, Academic Press, Inc. N.Y., (1990). 

Other suitable amplification methods include, but are not limited to 
polymerase chain reaction (PCR) (Innis, et aL, PCR Protocols. A guide to Methods 
and Application. Academic Press, Inc. San Diego, (1990)), ligase chain reaction 
(LCR) (see Wu and Wallace, Genomics, 4: 560 (1989), Landegren, et aL, Science, 
241: 1077 (1988) and Barringer, et aL, Gene, 89: 117 (1990), transcription 
ampliation (Kwoh,etal.,i^ C/&4, 86: 1173 (1989)), and self- 

sustained sequence replication (Guatelli, et aL, Proc. Nat Acad. Set USA, 87: 1874 
(1990)). 

Cell lysates or tissue homogenates often contain a number of inhibitors of 
polymerase activity. Therefore, RT-PCR typically incorporates preliminary steps 
to isolate total RNA or mRNA for subsequent use as an amplification template. 
A one-tube mRNA capture method may be used to prepare poly(A) + RNA samples 
suitable for immediate RT-PCR in the same tube (Boehringer Mannheim). The 
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captured mRNA can be directly subjected to RT-PCR by adding a reverse 
transcription mix and, subsequently, a PCR mix. 

In a particularly preferred embodiment, the sample mRNA is reverse 
transcribed with a reverse transcriptase and a primer consisting of ohgo(dT) and 
5 a sequence encoding the phage T7 promoter to provide single stranded DNA 

template. The second DNA strand is polymerized using a DNA polymerase. After 
synthesis of double-stranded cDNA, T7 RNA polymerase is added and RNA is 
transcribed from the cDNA template. Successive rounds of transcription from each 
single cDNA template results in amplified RNA. Methods of in vitro 

10 polymerization are well known to those of skill in the art (see, e.g. 9 Sambrook, 

supra.) and this particular method is described in detail by Van Gelder, et al, Proa 
Natl Acad. Sci. USA, 87: 1663-1667 (1990) who demonstrate that in vitro 
amplification according to this method preserves the relative frequencies of the 
various RNA transcripts. Moreover, Eberwine et al Proa Natl Acad Set USA, 

15 89: 3010-3014 provide a protocol that uses two rounds of amplification via in vitro 

transcription to achieve greater than 10^ fold amplification of the original starting 
material, thereby permitting expression monitoring even where biological samples 
are limited. 

It will be appreciated by one of skill in the art that the direct transcription 
20 method described above provides an antisense (aRNA) pool. Where antisense 

RNA is used as the target nucleic acid, the oligonucleotide probes provided in the 
array are chosen to be complementary to subsequences of the antisense nucleic 
acids. Conversely, where the target nucleic acid pool is a pool of sense nucleic 
acids, the oligonucleotide probes are selected to be complementary to subsequences 
25 of the sense nucleic acids. Finally, where the nucleic acid pool is double stranded, 

the probes may be of either sense as the target nucleic acids include both sense and 
antisense strands. 

The protocols cited above include methods of generating pools of either 
sense or antisense nucleic acids. Indeed, one approach can be used to generate 
30 either sense or antisense nucleic acids as desired. For example, the cDNA can be 

directionally cloned into a vector (e.g., Stratagene's p Bluscript II KS (+) 
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phagemid) such that it is flanked by the T3 and 17 promoters. In vitro 
transcription with the T3 polymerase will produce RNA of one sense (the sense 
depending on the orientation of the insert), while in vitro transcription with the 17 
polymerase will produce RNA having the opposite sense. Other suitable cloning 
systems include phage lambda vectors designed for Cre-/oxP plasmid subcloning 
{see e.g., Palazzolo etal, Gene, 88: 25-36 (1990)). 

(B) Hybridizing nucleic acids to high density arrays 
1. Probe design 

One of skill in the art will appreciate that an enormous number of array 
designs are suitable for the practice of this invention. The high density array will 
typically include a number of probes that specifically hybridize to the sequences 
of interest. In addition, in a preferred embodiment, the array will include one or 
more control probes. 

The high density array chip includes "test probes." Test probes could be 
oligonucleotides that range from about 5 to about 45 or 5 to about 500 nucleotides, 
more preferably from about 10 to about 40 nucleotides and most preferably from 
about 15 to about 40 nucleotides in length. In other particularly preferred 
embodiments the probes are 20 or 25 nucleotides in length. In another preferred 
embodiments, test probes are double or single strand DNA sequences. DNA 
sequences are isolated or cloned from nature sources or amplified from nature 
sources using nature nucleic acid as templates. These probes have sequences 
complementary to particular subsequences of the genes whose expression they are 
designed to detect Thus, the test probes are capable of specifically hybridizing to 
the target nucleic acid they are to detect. 

hi addition to test probes that bind the target nucleic acid(s) of interest, the 
high density array can contain a number of control probes. The control probes fall 
into three categories referred to herein as 1) normalization controls; 2) expression 
level controls; and 3) mismatch controls. 
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Normalization controls are oligonucleotide or other nucleic acid proves that 
are complementary to labeled reference oligonucleotides or other nucleic acid 
sequences that are added to the nucleic acid sample. The signals obtained from the 
normalization controls after hybridization provide a control for variations in 
5 hybridization conditions, label intensity, "reading" efficiency and other factors that 

may cause the signal of a perfect hybridization to vary between arrays. In a 
preferred embodiment, signals (e.g. y fluorescence intensity) read from all other 
probes in the array are divided by the signal (eg., fluorescence intensity) from the 
control probes thereby normalizing the measurements. 

10 Virtually any probe may serve as a normalization control. However, it is 

recognized that hybridization efficiency varies with base composition and probe 
length. Preferred normalization probes are selected to reflect the average length of 
the other probes present in the array, however, they can be selected to cover a range 
of lengths. The normalization controls) can also be selected to reflect the 

1 5 (average) base composition of the other probes in the array, however in a preferred 

embodiment, only one or a few normalization probes are used and they are selected 
such that they hybridize well (ie. no secondary structure) and do not match any 
target-specific probes. 

Expression level controls are probes that hybridize specifically with 

20 constitutively expressed genes in the biological sample. Virtually any 

constitutively expressed gene provides a suitable target for expression level 
controls. Typically expression level control probes have sequences complementary 
to subsequences of constitutively expressed "housekeeping genes" including, but 
not limited to the B-actin gene, the transferrin receptor gene, the GAPDH gene, and 

25 the like. 
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Mismatch controls may also be provided for the probes to the target genes, 
for expression level controls or for normalization controls. Mismatch controls are 
oligonucleotide probes or other nucleic acid probes identical to their corresponding 
test or control probes except for the presence of one or more mismatched bases. 
5 A mismatched base is a base selected so that it is not complementary to the 

corresponding base in the target sequence to which the probe would otherwise 
specifically hybridize. One or more mismatches are selected such that under 
appropriate hybridization conditions (e.g. stringent conditions) the test or control 
probe would be expected to hybridize with its target sequence, but the mismatch 

10 probe would not hybridize (or would hybridize to a significantly lesser extent). 

Preferred mismatch probes contain a central mismatch. Thus, for example, where 
a probe is a 20 mer, a corresponding mismatch probe will have the identical 
sequence except for a single base mismatch (e.g., substituting a G, a C or a T for 
an A) at any of positions 6 through 14 (the central mismatch). 

1 5 Mismatch probes thus provide a control for non-specific binding or cross- 

hybridization to a nucleic acid in the sample other than the target to which the 
probe is directed. Mismatch probes thus indicate whether a hybridization is specific 
or not. For example, if the target is present the perfect match probes should be 
consistently brighter than the mismatch probes. In addition, if all central 

20 mismatches are present, the mismatch probes can be used to detect a mutation. The 

difference in intensity between the perfect match and the mismatch probe (I(PM)- 
I(MM)) provides a good measure of the concentration of the hybridized material. 

The high density array may also include sample preparation/amplification 
control probes. These are probes that are complementary to subsequences of 

25 control genes selected because they do not normally occur in the nucleic acids of 

the particular biological sample being assayed. Suitable sample 
preparation/amplification control probes include, for example, probes to bacterial 
genes (e.g. 9 Bio B) where the sample in question is a biological from a eukaryote. 
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The RNA sample is then spiked with a known amount of the nucleic acid 
to which the sample preparation/amplification control probe is directed before 
processing. Quantification of the hybridization of the sample 
preparation/amplification control probe then provides a measure of alteration in the 
5 abundance of the nucleic acids caused by processing steps (e.g. PCR, reverse 

transcription, in vitro transcription, etc). 

In a preferred embodiment, oligonucleotide probes in the high density array 
are selected to bind specifically to the nucleic acid target to which they are directed 
with minimal non-specific binding or cross-hybridization under the particular 

1 0 hybridization conditions utilized. Because the high density arrays of this invention 

can contain in excess of 1,000,000 different probes, it is possible to provide every 
probe of a characteristic length that binds to a particular nucleic acid sequence. 
Thus, for example, the high density array can contain every possible 20-mer 
sequence complementary to an IL-2 mRNA. 

15 However, there may exist 20-mer subsequences that are not unique to the 

EL-2 mRNA. Probes directed to these subsequences are expected to cross- 
hybridize with occurrences of their complementary sequence in other regions of the 
sample genome. Similarly, other probes simply may not hybridize effectively 
under the hybridization conditions due to secondary structure, or interactions 

20 with the substrate or other probes). Thus, in a preferred embodiment, the probes 

that show such poor specificity or hybridization efficiency are identified and may 
not be included either in the high density array itself (eg., during fabrication of the 
array) or in the post-hybridization data analysis. 

In addition, in a preferred embodiment, expression monitoring arrays are 

25 used to identify the presence and expression (transcription) level of genes which 

are several hundred base pairs long. For most applications it would be useful to 
identify the presence, absence, or expression level of several thousand to one 
hundred thousand genes. Because the number of oligonucleotides per array is 
. limited in a preferred embodiment, it is desired to include only a limited set of 

30 probes specific to each gene whose expression is to be detected. 

23 



WO 00/011218 



PCT/US99/18772 



As disclosed in U.S. Application Ser. No. 08/772,376, probes as short as 
15, 20, or 25 nucleotide are sufficient to hybridize to a subsequence of a gene and 
that, for most genes, there is a set of probes that performs well across a wide range 
of target nucleic acid concentrations. In a preferred embodiment, it is desirable to 
5 choose a preferred or "optimum" subset of probes for each gene before 

synthesizing the high density array. 

2. Forming High Density Arrays. 

10 Methods of forming high density arrays of oligonucleotides, peptides and 

other polymer sequences with a minimal number of synthetic steps are known. The 
oligonucleotide analogue array can be synthesized on a solid substrate by a variety 
of methods, including, but not limited to, light-directed chemical coupling, and 
mechanically directed coupling. See Pirrung et al, U.S. Patent No. 5,143,854 (see 

15 also PCT Application No. WO 90/15070) and Fodor et aL, PCT Publication Nos. 

WO 92/10092 and WO 93/09668 and US Ser. No. 07/980,523 which disclose 
methods of forming vast arrays of peptides, oligonucleotides and other molecules 
using, for example, light-directed synthesis techniques. See also, Fodor et aL, 
Science, 251, 767-77 (1991). These procedures for synthesis of polymer arrays are 

20 now referred to as VLSIPS™ procedures. Using the VLSIPS™ approach, one 

heterogeneous array of polymers is converted, through simultaneous coupling at 
a number of reaction sites, into a different heterogeneous array. See, U.S. 
Application Serial Nos. 07/796,243 and 07/980,523. 
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The development of VLSIPS™ technology as described in the 
above-noted U.S. Patent No. 5,143,854 and PCT patent publication Nos. WO 
90/15070 and 92/10092, is considered pioneering technology in the fields of 
combinatorial synthesis and screening of combinatorial libraries. More recently, 
5 patent application Serial No. 08/082,937, filed June 25, 1993, describes methods 

for making arrays of oligonucleotide probes that can be used to check or determine 
a partial or complete sequence of a target nucleic acid and to detect the presence 
of a nucleic acid containing a specific oligonucleotide sequence. 

In brie£ the light-directed combinatorial synthesis of oligonucleotide arrays 

10 on a glass surface proceeds using automated phosphoramidite chemistry and chip 

masking techniques. In one specific implementation, a glass surface is derivatized 
with a silane reagent containing a functional group, eg., a hydroxyl or amine group 
blocked by a photolabile protecting group. Photolysis through a photolithogaphic 
mask is used selectively to expose functional groups which are then ready to react 

15 with incoming 5-photoprotected nucleoside phosphoramidites. The 

phosphoramidites react only with those sites which are illuminated (and thus 
exposed by removal of the photolabile blocking group). Thus, the 
phosphoramidites only add to those areas selectively exposed from the preceding 
step. These steps are repeated until the desired array of sequences have been 

20 synthesized on the solid surface. Combinatorial synthesis of different 

oligonucleotide analogues at different locations on the array is determined by the 
pattern of illumination during synthesis and the order of addition of coupling 
reagents. 

In the event that an oligonucleotide analogue with a polyamide backbone 
25 is used in the VLSIPS™ procedure, it is generally inappropriate to use 

phosphoramidite chemistry to perform the synthetic steps, since the monomers do 
not attach to one another via a phosphate linkage. Instead, peptide synthetic 
methods are substituted. See, e.g., Pirrung et aL U.S. Pat. No. 5,143,854. 

Peptide nucleic acids are commercially available from, e.g., Biosearch, Inc. 
30 (Bedford, MA) which comprise a polyamide backbone and the bases found in 

naturally occurring nucleosides. Peptide nucleic acids are capable of binding to 
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nucleic acids with high specificity, and are considered "oligonucleotide analogues" 
for purposes of this disclosure. 

In addition to the foregoing, additional methods which can be used to 
generate an array of oligonucleotides on a single substrate are described in co- 
5 pending Applications Ser. No. 07/980,523, filed November 20, 1992, and 

07/796,243, filed November 22, 1991 and in PCT Publication No. WO 93/09668. 
In the methods disclosed in these applications, reagents are delivered to the 
substrate by either (1) flowing within a channel defined on predefined regions or 
(2) "spotting 11 on predefined regions or (3) through the use of photoresist. 

10 However, other approaches, as well as combinations of spotting and flowing, may 

be employed. In each instance, certain activated regions of the substrate are 
mechanically separated from other regions when the monomer solutions are 
delivered to the various reaction sites. 

A typical "flow channel" method applied to the compounds and libraries of 

15 the present invention can generally be described as follows. Diverse polymer 

sequences are synthesized at selected regions of a substrate or solid support by 
forming flow channels on a surface of the substrate through which appropriate 
reagents flow or in which appropriate reagents are placed. For example, assume 
a monomer "A" is to be bound to the substrate in a first group of selected regions. 

20 If necessary, all or part of the surface of the substrate in all or a part of the selected 

regions is activated for binding by, for example, flowing appropriate reagents 
through all or some of the channels, or by washing the entire substrate with 
appropriate reagents. After placement of a channel block on the surface of the 
substrate, a reagent having the monomer A flows through or is placed in all or 

25 some of the channels). The channels provide fluid contact to the first selected 

regions, thereby binding the monomer A on the substrate directly or indirectly (via 
a spacer) in the first selected regions. 

Thereafter, a monomer B is coupled to second selected regions, some of 
which may be included among the first selected regions. The second selected 

30 regions will be in fluid contact with a second flow channel(s) through translation, 

rotation, or replacement of the channel block on the surface of the substrate; 
through opening or closing a selected valve; or through deposition of a layer of 
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chemical or photoresist If necessary, a step is performed for activating at least the 
second regions. Thereafter, the monomer B is flowed through or placed in the 
second flow channel(s), binding monomer B at the second selected locations. In 
this particular example, the resulting sequences bound to the substrate at this stage 
of processing will be, for example, A, B, and AB. The process is repeated to form 
a vast array of sequences of desired length at known locations on the substrate. 

After the substrate is activated, monomer A can be flowed through some 
of the channels, monomer B can be flowed through other channels, a monomer C 
can be flowed through still other channels, etc. In this manner, many or all of the 
reaction regions are reacted with a monomer before the channel block must be 
moved or the substrate must be washed and/or reactivated. By making use of many 
or all of the available reaction regions simultaneously, the number of washing and 
activation steps can be minimized. 

One of skill in the art will recognize that there are alternative methods of 
framing channels or otherwise protecting a portion of the surface of the substrate. 
For example, according to some embodiments, a protective coating such as a 
hydrophilic or hydrophobic coating (depending upon the nature of the solvent) is 
utilized over portions of the substrate to be protected, sometimes in combination 
with materials that facilitate wetting by the reactant solution in other regions. In 
this manner, the flowing solutions are further prevented from passing outside of 
their designated flow paths. 

High density nucleic acid arrays can be fabricated by depositing 
presynthesized or natural nucleic acids in predetermined positions. Synthesized 
or natural nucleic acids are deposited on specific locations of a substrate by light 
directed targeting and oligonucleotide directed targeting. Nucleic acids can also 
be directed to specific locations in much the same manner as the flow channel 
methods. For example, a nucleic acid A can be delivered to and coupled with a 
first group of reaction regions which have been appropriately activated. Thereafter, 
a nucleic acid B can be delivered to and reacted with a second group of activated 
reaction regions. Nucleic acids are deposited in selected regions. Another 
embodiment uses a dispenser that moves from region to region to deposit nucleic 
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acids in specific spots. Typical dispensers include a micropipette or capillary pin 
to deliver nucleic acid to the substrate and a robotic system to control the position 
of the micropipette with respect to the substrate. In other embodiments, the 
dispenser includes a series of tubes, a manifold, an array of pipettes or capillary 
5 pins, or the like so that various reagents can be delivered to the reaction regions 

simultaneously. 

3. Hybridization 

Nucleic acid hybridization simply involves contacting a probe and target 

10 nucleic acid under conditions where the probe and its complementary target can 

form stable hybrid duplexes through complementary base pairing. The nucleic 
acids that do not form hybrid duplexes are then washed away leaving the 
hybridized nucleic acids to be detected, typically through detection of an attached 
detectable label. It is generally recognized that nucleic acids are denatured by 

15 increasing the temperature or decreasing the salt concentration of the buffer 

containing the nucleic acids. Under low stringency conditions low 
temperature and/or high salt) hybrid duplexes (e.g 9 DNArDNA, RNA:RNA, or 
RNA:DNA) will form even where the annealed sequences are not perfectly 
complementary. Thus specificity of hybridization is reduced at lower stringency. 

20 Conversely, at higher stringency (e.g., higher temperature or lower salt) successful 

hybridization requires fewer mismatches. 

One of skill in the art will appreciate that hybridization conditions may be 
selected to provide any degree of stringency. In a preferred embodiment, 
hybridization is performed at low stringency in this case in 6X SSPE-T at 37 C 

25 (0.005% Triton X-100) to ensure hybridization and then subsequent washes are 

performed at higher stringency (e.g , 1 X SSPE-T at 37 C) to eliminate mismatched 
hybrid duplexes. Successive washes may be performed at increasingly higher 
stringency (e.g 9 down to as low as 0.25 X SSPE-T at 37 C to 50 C) until a desired 
level of hybridization specificity is obtained. Stringency can also be increased by 

30 addition of agents such as fonnamide. Hybridization specificity may be evaluated 

by comparison of hybridization to the test probes with hybridization to the various 
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controls that can be present (e.g, expression level control, normalization control, 

mismatch controls, etc.). 

In general, there is a tradeoff between hybridization specificity (stringency) 

and signal intensity. Thus, in a preferred embodiment, the wash is performed at the 
5 highest stringency that produces consistent results and that provides a signal 

intensity greater than approximately 10% of the background intensity. Thus, in a 

preferred embodiment, the hybridized array may be washed at successively higher 

stringency solutions and read between each wash. Analysis of the data sets thus 

produced will reveal a wash stringency above which the hybridization pattern is not 
10 appreciably altered and which provides adequate signal for the particular 

oligonucleotide probes of interest. 

In a preferred embodiment, background signal is reduced by the use of a 

detergent (e.g. 9 C-TAB) or a blocking reagent (e.g 9 sperm DNA, cot-1 DNA, etc.) 

during the hybridization to reduce non-specific binding. In a particularly preferred 
15 embodiment, the hybridization is performed in the presence of about 0.5 mg/ml 

DNA (e.g. 9 herring sperm DNA). The use of blocking agents in hybridization is 

well known to those of skill in the art (see, e.g 9 Chapter 8 in P. Tijssen, supra.) 
The stability of duplexes formed between RNAs or DNAs are generally in 

the order of RNA:RNA > KNA:DNA > DNA:DNA, in solution. Long probes have 
20 better duplex stability with a target, but poorer mismatch discrimination than 

shorter probes (mismatch discrimination refers to the measured hybridization 

signal ratio between a perfect match probe and a single base mismatch probe). 

Shorter probes (e.g., 8-mers) discriminate mismatches very well, but the overall 

duplex stability is low. 

25 Altering the thermal stability (1^ of the duplex formed between the target 

and the probe using, e.g. f known oligonucleotide analogues allows for optimization 
of duplex stability and mismatch discrimination. One useful aspect of altering the 
T m arises from the fact that adenine-thymine (A-T) duplexes have a lower T m than 
guanine-cytosine (G-C) duplexes, due in part to the fact that the A-T duplexes have 

30 2 hydrogen bonds per base-pair, while the G-C duplexes have 3 hydrogen bonds 

per base pair. In heterogeneous oligonucleotide arrays in which there is a non- 
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uniform distribution of bases, it is not generally possible to optimize hybridization 
for each oligonucleotide probe simultaneously. Thus, in some embodiments, it is 
desirable to selectively destabilize G-C duplexes and/or to increase the stability of 
A-T duplexes. This can be accomplished, e.g, by substituting guanine residues in 
5 the probes of an array which form G-C duplexes with hypoxanthine, or by 

substituting adenine residues in probes which form A-T duplexes with 2,6 
diaminopurine or by using the salt tetramethyl ammonium chloride (TMAC1) in 
place of NaCl. 

Altered duplex stability conferred by using oligonucleotide analogue probes 

10 can be ascertained by following, eg., fluorescence signal intensity of 

oligonucleotide analogue arrays hybridized with a target oligonucleotide over time. 
The data allow optimization of specific hybridization conditions at, e.g., room 
temperature (for simplified diagnostic applications in the future). 

Another way of verifying altered duplex stability is by following the signal 

15 intensity generated upon hybridization with time. Previous experiments using 

DNA targets and DNA chips have shown that signal intensity increases with time, 
and that the more stable duplexes generate higher signal intensities faster than less 
stable duplexes. The signals reach a plateau or "saturate" after a certain amount of 
time due to all of the binding sites becoming occupied. These data allow for 

20 optimization of hybridization, and determination of the best conditions at a 

specified temperature. 

Methods of optimizing hybridization conditions are well known to those of 
skill in the art (see, e.g. t Laboratory Techniques in Biochemistry and Molecular 
Biology, Vol 24: Hybridization With Nucleic Acid Probes \ P. Tijssen, ed. Elsevier, 

25 N.Y.,(1993)). 

(C) Signal Detection 

In a preferred embodiment, the hybridized nucleic acids are detected by 
detecting one or more labels attached to the sample nucleic acids. The labels may 
be incorporated by any of a number of means well known to those of skill in the 
30 art. However, in a preferred embodiment, the label is simultaneously incorporated 

30 
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during the amplification step in the preparation of the sample nucleic acids. Thus, 
for example, polymerase chain reaction (PCR) with labeled primers or labeled 
nucleotides will provide a labeled amplification product. In a preferred 
embodiment, transcription amplification, as described above, using a labeled 
5 nucleotide (e.g. fluorescein-labeled U1P and/or CTP) incorporates a label into the 

transcribed nucleic acids. 

Alternatively, a label may be added directly to the original nucleic acid 
sample (e.g., mRNA, polyA mRNA, cDNA, etc.) or to the amplification product 
after the amplification is completed Means of attaching labels to nucleic acids are 

10 well known to those of skill in the art and include, for example nick translation or 

end-labeling (e.g. with a labeled RNA) by kinasing of the nucleic acid and 
subsequent attachment (ligation) of a nucleic acid linker joining the sample nucleic 
acid to a label (e.g. , a fluorophore). 

Detectable labels suitable for use in the present invention include any 

15 composition detectable by spectroscopic, photochemical, biochemical, 

immunochemical, electrical, optical or chemical means. Useful labels in the 
present invention include biotin for staining with labeled streptavidin conjugate, 
magnetic beads (e.g, Dynabeads™), fluorescent dyes (e.g., fluorescein, texas red, 
rhodamine, green fluorescent protein, and the like), radiolabels (e.g., 3 H , 12 5 I, 

20 -*i>S, I4Q or ^2P\ enzymes (e.g., horse radish peroxidase, alkaline phosphatase 

and others commonly used in an ELIS A), and colorimetric labels such as colloidal 
gold or colored glass or plastic (e.g. 9 polystyrene, polypropylene, latex, etc.) beads. 
Patents teaching the use of such labels include U.S. Patent Nos. 3,817,837; 
3,850,752; 3,939,350; 3,996,345; 4,277,437; 4,275,149; and 4,366,241. 

25 Means of detecting such labels are well known to those of skill in the art. 

Thus, for example, radiolabels may be detected using photographic film or 
scintillation counters, fluorescent markers may be detected using a photodetector 
to detect emitted light Enzymatic labels are typically detected by providing the 
enzyme with a substrate and detecting the reaction product produced by the action 

30 of the enzyme on the substrate, and colorimetric labels are detected by simply 
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visualizing the colored label. One particular preferred methods uses colloidal gold 
label that can be detected by measuring scattered light. 

The label may be added to the target (sample) nucleic acid(s) prior to, or 
after the hybridization. So called "direct labels" are detectable labels that are 
directly attached to or incorporated into the target (sample) nucleic acid prior to 
hybridization. In contrast, so called "indirect labels" are joined to the hybrid 
duplex after hybridization. Often, the indirect label is attached to a binding moiety 
that has been attached to the target nucleic acid prior to the hybridization. Thus, 
for example, the target nucleic acid may be biotinylated before the hybridization. 
After hybridization, an avidin-conjugated fluorophore will bind the biotin bearing 
hybrid duplexes providing a label that is easily detected. For a detailed review of 
methods of labeling nucleic acids and detecting labeled hybridized nucleic acids 
see Laboratory Tecltniques in Biochemistry and Molecular Biology, Vol 24: 
Hybridization With Nucleic Acid Probes, P. Tijssen, ed. Elsevier, N.Y., (1993)). 

Fluorescent labels are preferred and easily added during an in vitro 
transcription reaction. In a preferred embodiment, fluorescein labeled UTP and 
CTP are incorporated into the RNA produced in an in vitro transcription reaction 
as described above. 

Means of detecting labeled target (sample) nucleic acids hybridized to the 
probes of the high density array are known to those of skill in the art. Thus, for 
example, where a colorimetric label is used, simple visualization of the label is 
sufficient. Where a radioactive labeled probe is used, detection of the radiation 
(e.g. with photographic film or a solid state detector) is sufficient. 

In a preferred embodiment, however, the target nucleic acids are labeled 
with a fluorescent label and the localization of the label on the probe array is 
accomplished with fluorescent microscopy. The hybridized array is excited with 
a light source at the excitation wavelength of the particular fluorescent label and 
the resulting fluorescence at the emission wavelength is detected. In a particularly 
preferred embodiment, the excitation light source is a laser appropriate for the 
excitation of the fluorescent label. 

The confocal microscope may be automated with a computer-controlled 
stage to automatically scan the entire high density array. Similarly, the microscope 
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may be equipped with a phototransducer (e.g., a photomultiplier, a solid state array, 
a CCD camera, etc.) attached to an automated data acquisition system to 
automatically record the fluorescence signal produced by hybridization to each 
oligonucleotide probe on the array. Such automated systems are described at 
5 length in U.S. Patent No: 5,143,854, PCT Application 20 92/10092, and copending 

U.S. Application Ser. No. 08/195,889 filed on February 10, 1994. Use of laser 
illumination in conjunction with automated confocal microscopy for signal 
detection permits detection at a resolution of better than about 100 ^m, more 
preferably better than about 50 jam, and most preferably better than about 25 ^m. 

10 One of skill in the art will appreciate that methods for evaluating the 

hybridization results vary with the nature of the specific probe nucleic acids used 
as well as the controls provided. la the simplest embodiment, simple 
quantification of the fluorescence intensity for each probe is determined. This is 
accomplished simply by measuring probe signal strength at each location 

15 (representing a different probe) on the high density array (e.g., where the label is 

a fluorescent label, detection of the amount of florescence (intensity) produced by 
a fixed excitation illumination at each location on the array). Comparison of the 
absolute intensities of an array hybridized to nucleic acids from a "test" sample 
with intensities produced by a "control" sample provides a measure of the relative 

20 expression of the nucleic acids that hybridize to each of the probes. 

One of skill in the art, however, will appreciate that hybridization signals 
will vary in strength with efficiency of hybridization, the amount of label on the 
sample nucleic acid and the amount of the particular nucleic acid in the sample. 
Typically nucleic acids present at very low levels (e.g. 9 < lpM) will show a very 

25 weak signal. At some low level of concentration, the signal becomes virtually 

indistinguishable from background. In evaluating the hybridization data, a 
threshold intensity value may be selected below which a signal is not counted as 
being essentially indistinguishable from background. 

Where it is desirable to detect nucleic acids expressed at lower levels, a 

30 lower threshold is chosen. Conversely, where only high expression levels are to 
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be evaluated a higher threshold level is selected. Ia a preferred embodiment, a 
suitable threshold is about 10% above that of the average background signal. 

In addition, the provision of appropriate controls permits a more detailed 
analysis that controls for variations in hybridization conditions, cell health, non- 
specific binding and the like. Thus, for example, in a preferred embodiment, the 
hybridization array is provided with normalization controls. These normalization 
controls are probes complementary to control sequences added in a known 
concentration to the sample. Where the overall hybridization conditions are poor, 
the normalization controls will show a smaller signal reflecting reduced 
hybridization. Conversely, where hybridization conditions are good, the 
normalization controls will provide a higher signal reflecting the improved 
hybridization. Normalization of the signal derived from other probes in the array 
to the normalization controls thus provides a control for variations in hybridization 
conditions. Typically, normalization is accomplished by dividing the measured 
signal from the other probes in the array by the average signal produced by the 
normalization controls. Normalization may also include correction for variations 
due to sample preparation and amplification. Such normalization may be 
accomplished by dividing the measured signal by the average signal from the 
sample preparation/amplification control probes {e.g., the Bio B probes). The 
resulting values may be multiplied by a constant value to scale the results. 

As indicated above, the high density array can include mismatch controls. 
In a preferred embodiment, there is a mismatch control having a central mismatch 
for every probe (except the normalization controls) in the array. It is expected that 
after washing in stringent conditions, where a perfect match would be expected to 
hybridize to the probe, but not to the mismatch, the signal from the mismatch 
controls should only reflect non-specific binding or the presence in the sample of 
a nucleic acid that hybridizes with the mismatch. Where both the probe in question 
and its corresponding mismatch control both show high signals, or the mismatch 
shows a higher signal than its corresponding test probe, there is a problem with the 
hybridization and the signal from those probes is ignored. The difference in 
hybridization signal intensity between the target specific probe and its 
corresponding mismatch control is a measure of the (ftscrimination of the target- 
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specific probe. Thus, in a preferred embodiment, the signal of the mismatch probe 
is subtracted from the signal from its corresponding test probe to provide a measure 
of the signal due to specific binding of the test probe. 

The concentration of a particular sequence can then be determined by 
measuring the signal intensity of each of the probes that bind specifically to that 
gene and normalizing to the normalization controls. Where the signal from the 
probes is greater than the mismatch, the mismatch is subtracted. Where the 
mismatch intensity is equal to or greater than its corresponding test probe, the 
signal is ignored. The expression level of a particular gene can then be scored by 
the number of positive signals (either absolute or above a threshold value), the 
intensity of the positive signals (either absolute or above a selected threshold 
value), or a combination of both metrics (e.g. 9 a weighted average). 

In some preferred embodiments, a computer system is used to compare the 
hybridization intensities of the perfect match and mismatch probes of each pair. 
If the gene is expressed, the hybridization intensity (or affinity) of a perfect match 
probe of a pair should be recognizably higher than the corresponding mismatch 
probe. Generally, if the hybridizations intensities of a pair of probes are 
substantially the same, it may indicate the gene is not expressed. However, the 
detenmnation is not based on a single pair of probes, the determination of whether 
a gene is expressed is based on an analysis of many pairs of probes. 

After the system compares the hybridization intensity of the perfect 
match and mismatch probes, the system indicates expression of the gene. As an 
example, the system may indicate to a user that the gene is either present 
(expressed), marginal or absent (unexpressed). Specific procedures for data 
analysis is disclosed in U.S. Application 08/772,376, previously incorporated 
for all purposes. 

In addition to high density nucleic acid arrays, other methods are also 
useful for massive gene expression monitoring. Differential display, described 
by Liang, P. and Pardee, A.B. (Differential Display of eukaryotic messenger 
RNA by means of the polymerase chain reaction. Science 257:967-971, 1992, 
incorporated herein by reference for all purposes) provides a useful mean for 
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distinguishing gene expression between two samples. Serial analysis of gene 
expression, described by Velculescu et al. (Serial Analysis of Gene Expression. 
Science, 270:484-487, 1995, incorporated herein by reference for all purposes) 
provides another method for quantitative and qualitative analysis of gene 
expression. Optical fiber oligonucleotide sensors, described by Ferguson et al. 
(A Fiber-optic DNA biosensor microarray for the analysis of gene expression. 
Nature-Biotechnology 14:1681-1684, 1996), can also be used for gene 
expression monitoring. 

It is understood that the examples and embodiments described herein are 
for illustrative purposes only and that various modifications or changes in light 
thereof will be suggested to persons skilled in the art and are to be included 
within the spirit and purview of this application and scope of the appended 
claims. All publications, patents, and patent applications cited herein are hereby 
incorporated by reference for all purposes. 

EXAMPLES 

SUMMARY 

Mechanistic insights to viral replication and pathogenesis generally 
have come from the analysis of viral gene products, either by studying their 
biochemical activities and interactions individually or by creating mutant 
viruses and analyzing their phenotype. Now it is possible to identify and 
catalog the host cell genes whose mRNA levels change in response to the 
pathogen* We have employed DNA array technology to monitor the level of 
approximately 6600 human mRNAs in uninfected as compared to human 
cytomegalovirus-infected cells. The level of 258 mRNAs changed by a 
factor of 4 or more early after infection, before the onset of viral DNA 
replication. Several of these mRNAs encode gene products that might play 
key roles in virus-induced pathogenesis, identifying them as intriguing 
targets for further study. 
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INTRODUCTION 

Human cytomegalovirus (HCMV) has the potential to alter cellular gene 
expression though multiple mechanisms. Its initial interaction with the cell 
surface could initiate a regulatory signal; indeed, the virion gB and gH 
glycoproteins induce cellular transcription factors when added to uninfected 
cells (1). Constituents of the virion, such as the tegument protein, pp71, migrate 
to the nucleus and activate transcription after infection (2), and viral proteins 
synthesized after infection, such as the immediate early 1 and 2 proteins, 
modulate transcription (3-5). The virus encodes several G protein-coupled 
receptors (6,7) that likely initiate gene regulatory signal cascades in response to 
ligands, and HCMV infection has been shown to perturb cell cycle regulation 
(8-11), which leads to changes in cellular gene expression. The complex 
virus-host cell interaction has the potential to dramatically modulate the 
expression of cellular genes. 

Relatively few cellular genes have been identified whose activity 
changes in HCMV-infected cells (12). Recently, we used differential display 
analysis to identify 15 interferon-inducible genes that are activated by the virus 
subsequent to infection (13). However, this screen identified only genes whose 
mRNA levels changed dramatically, and we were not able to perform the screen 
under a variety of conditions given its labor-intensive nature. In contrast to 
differential display, the DNA array assay is easily performed and can detect 
subtle changes in mRNA levels. We report the identification of 258 cellular 
mRNAs whose level changes by a factor of 4 or more before the onset of 
HCMV DNA replication. 

MATERIALS AND METHODS 
Cells and Viruses. Primary human foreskin fibroblasts at passage 10-15 were 
cultured in DMEM containing 10% fetal calf serum. After the cells remained at 
confluence for thee days, they were infected at a multiplicity of 3 
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plaque-forming units per cell with HCMV AD169 or Toledo virions that were 
purified as described (14). 

Sample preparation and analysis with DNA arrays. Biotinylated 
single-stranded antisense RNA samples for hybridization were prepared as 
5 described (1 5) with minor modifications. Total cellular RNA was prepared 

using the TRIZOL Reagent (GibcoBRL), polyadenylated RNA was isolated, 
and portions (5 jag) were used as the template for the first strand cDNA 
synthesis in a reaction that was primed with oligo (dT) containing a T7 RNA 
polymerase promoter sequence at its 5' end 

10 [5 , -GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG(T) 24 -3 , 3. 

The second cDNA strand was synthesized using E. coli DNA polymerase I and 
ligase. The resulting cDNA (0.5-1 fig) was used as template to make a 
biotinylated RNA probe by in vitro transcription using the T7 Megascript 
System (Ambion). Unincorporated nucleotides were removed using a G-50 

15 Quick Spin Column (Boehringer Mannheim). The labeled RNA was fragmented 

to an average size of 50-100 bases by incubating at 94°C for 30 min in buffer 
containing 40 mM Tris-Ac, pH8.1, 100 mM KOAc, and 30 mM MgOAc. The 
hybridization (15 h), washing and staining protocols were as described (15), and 
employed a set of four human gene chips (HUM6000 A, B, C and D, 

20 Asymetrix). The DNA arrays were scanned using a confocal scanner 

manufactured for Afiymetrix by Molecular Dynamics. 

Data analysis: The data collected in each hybridization experiment was 
processed using the GeneChip ™ software supplied with the Afiymetrix 
instrumentation system. To evaluate if RNA corresponding to each of the 6600 

25 genes encoded on the array was detectable or undetectable a number of 

parameters were evaluated (15, 16), including the number of probe pairs 
interrogating each gene in which the intensity of the perfect match (PM)- 
hybridization reaction exceeded that of the mismatch (MM) hybridization signal 
(cutoffs <45%) and the PM/MM ratios for each set of probe pairs. To determine 

30 the quantitative amounts of RNA from each gene, the average of the differences 

(PM-MM) for each probe pair in a probe set was calculated (cutoff >50) as well 
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as the average differences across the prohe sets (cutoff >10). The cutoff 
thresholds were determined empirically to be conservative, i.e., they minimized 
false positives. The change in the level of expression for any gene was 
considered significant if the change in the average differences across the probe 
5 sets was greater than 3-fold. 

RNA analysis by Northern blot GeneChip™ results were confirmed by 
Northern blot assay. Total RNA (3 fig) from mock-infected cells or cells 
infected with the HCMV AD169 or Toledo strains was subjected to 
electrophoresis, blotted to a membrane and probed with random 
1 0 hexanucleotide-primed 32 P-labeled cDNA fragments from I.M. A.G.E. clones 

(Genome Systems, Inc.). 
RESULTS 

The gene chip assay utilized a set of four probe arrays that together 
include oligonucleotides corresponding to more than 6600 human mRNAs (16). 

15 Each array (1 .6 cm 2 ) contains more than 65,000 features, and a different 

oligodeoxyribonucleotide (25 bases) is synthesized on the surface of the 
derivatized glass wafer within the boundaries of each feature using light 
sensitive chemistry (17-21). The arrays contain 20 pairs of oligonucleotide 
probes corresponding to each RNA that is interrogated. Each probe pair consists 

20 of one 25-mer that is a perfect complement to the RNA (a perfect match probe) 

and a companion oligonucleotide that carries a single base difference in a 
central position (a mismatch probe). The mismatch probes serve as internal 
controls for hybridization specificity. Empirically derived rules used for the 
selection of oligonucleotide probes with the best sensitivity and specificity have 

25 been described (16). 

RNA samples were prepared for analysis at 40 min, 8 h and 24 h after 
mock-infection or HCMV (strain AD169) infection of primary human 
fibroblasts. Under these conditions, HCMV DNA replication begins between 
24-36 h after infection (10), and the complete viral replication cycle requires 
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mRNA differed by more than a factor of 3 in the replicate experiments (Fig. 
IB). This control demonstrates that the hybridization signals observed in 
independent experiments are highly reproducible. Further, the two preparations 
of infected cell target RNAs were prepared from infected primary fibroblasts 
derived from two different tissue samples, ruling out the possibility that changes 
in RNA levels might reflect genetic differences in the host cells. Differences 
greater than 3-fold observed for hybridization signals in comparisons of 
mock-infected versus infected cells should identify genes whose mRNA levels 
change after infection. 

When target RNA preparations were compared at 40 min after mock or 
virus infection, the level of 27 mRNAs had changed in response to infection by 
a factor of thee or more; at 8 and 24 h after infection, the number of altered 
mRNAs increased to 93 and 364, respectively (Fig. 2). Applying a more 
stringent four-fold cut off, we generated a set of 258 mRNAs for further 
analysis (Table 1). Of these mRNAs, 124 increased and 134 decreased after 
infection. We assume that most changes result from altered transcriptional 
regulation, but we have not yet tested this supposition. We confirmed 49 (40%) 
of the mRNAs predicted to be increased and 23 (17%) of the mRNAs predicted 
to be decreased either by northern blot analysis of independent RNA 
preparations (representative results in Fig. 3) or by reference to earlier studies 
that demonstrated a change. All attempts to confirm a predicted alteration in the 
group of 258 mRNAs were successful. 

We assayed changes in mRNA levels for a total of 58 genes in this study 
by northern blot. When we performed these assays, we included RNA 
preparations from cells infected with HCMV strain AD 169, the laboratory 
adapted strain used for the DNA array analysis, and HCMV strain Toledo, a 
clinical isolate that has not been extensively passaged in cultured cells (22). We 
observed the same alteration in mRNA level for both infections (representative 
results in Fig. 3). Although we might find some differences as more genes are 
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HCMV alter cellular gene expression in a similar fashion. 
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DISCUSSION 

HCMV replicates in many different cell types within its infected host, 
some of which might respond to infection differently than the primary 
fibroblasts we have studied here. Keeping this caveat in mind, we nevertheless 
can speculate that several of the cellular genes whose mRNA levels change after 
infection of fibroblasts might profoundly influence HCMV replication and 
pathogenesis. 

HLA-E mRNAs. In order to protect infected cells from cytotoxic T 
lymphocytes, multiple HCMV gene products act to reduce cell surface 
expression of classical class I MHC molecules (23-28). Although these viral 
activities protect infected cells from cytotoxic T lymphocytes, they also have 
the potential to render infected cells susceptible to natural killer (NK) cells that 
can recognize and destroy cells that no longer express class I MHC molecules. 
HLA-E mRNA is induced by a factor of 19 at 24 h after infection (Table 1), 
whereas HLA-A, HLA-D and HLA-G family members that were represented in 
the DNA arrays were not changed (data not shown). HLA-E is a nonclassical 
class I molecule whose cell surface expression requires that it bind peptides 
derived from the signal sequences of other class I molecules (HLA-A, -B and 
-C) (29). Recently, it has been shown that natural killer (NK) cells recognize 
and spare target cells expressing HLA-E on their surface (30, 31). This 
recognition is mediated by the NK cell CD94-NKG2 cell surface receptor. 
Assuming that the elevated mRNA leads to elevated cell surface expression of 
HLA-E, this modulation should protect virus-infected cells from NK cell 
killing. This would be the second mechanism by which HCMV avoids NK cell 
surveillance. The viral UL18 protein is an MHC homologue that engages 
another receptor (NKIR) on the NK cell to avoid attack (32). 

Ro/SSA 52 kDa mRNA. HCMV-infected cells contain enhanced levels 
of the Ro/SSA 52 kDa protein (Table 1). This protein, which is a constituent of 
a ribonucleoprotein complex, is induced by a factor of 12 at 24 h after infection. 
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Autoantibodies to this protein are found in a variety of connective tissue 
diseases: commonly in systemic lupus erythematosis, neonatal lupus 
erythematosis, and Sjogren's syndrome, and less frequently in rheumatoid 
arthitis (33). There is good evidence that these autoantibodies play a direct 
pathogenic role in neonatal lupus erythematosis and subacute cutaneous lupus 
erythematosis (33, 34). However, the mechanism by which the immune system 
initially responds to Ro/JSSA and other intracellular self-antigens is not clear. 
One popular hypothesis suggests that molecular mimicry is an important 
initiating mechanism, i.e., aspects of the immune response to a microbe cross 
react with self-proteins (35). Conceivably, overexpression of a commonly 
targeted autoantigen, such as the Ro/SSA antigen in HCMV-infected cells, also 
could favor an autoimmune response. Although the Ro/SSA 52kDa antigen is 
normally found in the nucleus and cytoplasm, it can be detected on the surface 
of peripheral lymphocytes that have been stressed by heat shock or treatment 
with ultraviolet light (36). Perhaps stress induced by HCMV infection also leads 
to cell surface presentation of Ro/SSA, facilitating an autoimmune response to 
the overexpressed antigen. Murine cytomegalovirus has been shown to induce 
autoimmune antibodies in infected mice (37-40), although Ro/SSA antibodies 
were not monitored in these studies. 

Lipocortin 1, cPLA2 and COX-2 mRNAs. Multiple constituents of the 
pathway that produces prostaglandin E2 from arachidonic acid are modulated by 
HCMV (Table 1). Cytosolic phospholipase A2 (cPLA2) mRNA increases by a 
factor of 12 and cyclooxygenase-2 (COX-2) mRNA is elevated by a factor of 7 
at 24 h after infection. Lipocortin 1, also known as annexin I, mRNA decreases 
by a factor of 9 at 24 h after infection. When cPLA2 is activated by 
phosphorylation, it translocates to membranes where it selectively cleaves and 
releases arachidonic acid; then COX2 converts it to prostaglandin E2. 
Lipocortin 1 inhibits the activation of cPLA2 (41). Thus, in HCMV-infected 
fibroblasts, the synthesis of prostaglandin E2 is activated by the induction of 
cPLA2 and COX2 and the inhibition of the negative regulator lipocortin 1, 
assuming that the changes in mRNA levels translate to changes in active 
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proteins. Further, HCMV infection has been shown to activate latent cPLA2 by 
inducing its phosphorylation (42). Thus, this pathway is strongly induced at 
both the transcriptional and posttranslational levels after infection, and this 
should lead to a marked increase in the production of prostaglandin E2. 
Prostaglandins serve as second messengers to stimulate a variety of responses, 
including inflammation. Perhaps the activation of this pathway is a cellular 
reaction to HCMV infection designed to induce a cell-mediated response that 
will kill the infected cell and thereby inhibit spread of the infection. 
Alternatively, one might speculate that the virus either induces the pathway or 
fails to antagonize the induction as a strategy to facilitate spread of the virus 
within the infected host. Inflammation might serve to lure monocytes and 
monocytic precursors to the vicinity of the infected cells where they can be 
infected. Cells of the monocytic lineage harbor HCMV on a long-term basis in a 
latent state (43-45). 

It is possible that the concerted changes in cPLA2, COX-2 and 
lipocortin 1 are an indirect effect of HCMV gene action. IL-ip has been shown 
to regulate this set of genes (46) in the same manner as seen in infected cells. 
Although several reports have suggested that IL-1 p activity is decreased in 
cultures of HCMV-infected monocytes (47, 48), the HCMV EE1 gene has been 
shown to induce the accumulation of IL-1 p mRNA in transfected monocytes 
(49, 50). The IL-ip gene was not included in the oligonucleotide arrays assayed 
in this report, so we do not know if its mRNA is induced by infection of 
fibroblasts. 

Thombospondin-1 mRNA. Thombospondin-1 is a calcium-binding protein 
released upon platelet activation (51). It is a constituent of the extracellular 
matrix that regulates cell growth and differentiation, and it might potentiate 
tumor progression (52). Recently, thombospondin-l-deficient mice have been 
produced (53) whose lungs exhibit acute and chonic cell infiltrates with 
increased fibroblastic and epithelial cell proliferation, matrix deposition and 
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diffuse alveolar hemorrhage characteristic of pneumonia. HCMV causes a 
21-fold reduction in the level of thombospondin-1 mRNA by 24 h after 
infection (Table 1). Replication in the lung that leads to pneumonia is one of the 
principle consequences of active HCMV infection in immunosuppressed 
5 individuals (54). Given the phenotype of thombospondin 1 -deficient mice, one 

can speculate that the reduction in this mRNA might contribute to pneumonia 
induced by acute HCMV infection. 

MTTF mRNA. The microphthalmia-associated transcription factor (MTTF) is 
the product of the microphthalmia gene. Mice have been described with a 

10 variety of mutations in this gene (55), and the most severe manifestations of the 

mutations include microphthalmia, oeteopetrosis and dea&ess. In the human, 
MTTF mutations were identified in two families afflicted with Waardenberg 
syndrome type 2, which causes hearing loss and patchy pigmentation of the 
eyes, hair and skin (56). Infection of humans with HCMV early in pregnancy 

1 5 has been reported to cause anophthalmia (57) and congenital infection of mice 

with murine cytomegalovirus can cause microphthalmia (58). Modulation of 
MTTF mRNA levels by the virus could contribute to these abnormalities. MTTF 
mRNA is reduced by a factor of 4-8 at 24 h after HCMV infection of 
fibroblasts. Whereas the association of HCMV with eye abnormalities appears 

20 to be rare, congenital HCMV infection is a common cause of hearing loss. 

Conceivably, HCMV-induced hearing loss is a consequence of an inhibitory 
effect on MTTF mRNA expression during development. This supposition is 
consistent with the observation that MTTF mutations are associated with hearing 
loss in the Waardenberg syndrome. HCMV could potentially modulate MTTF in 

25 cells that are eventually killed or in cells where viral gene expression does not 

lead to cell death. 

Conclusion. The roles of the cellular genes discussed above in HCMV 
replication and pathogenesis remain highly speculative. Nevertheless, the ability 
to identify cellular genes whose functions provide tantalizing hints of potential 
30 mechanistic roles in infectious disease processes underscores the utility of gene 

array technology in the study of pathogens. The global analysis of changes in 
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mRNA levels provides a catalog of genes that are modulated as a result of the 
host-pathogen interaction and therefore deserve further scrutiny. DNA array 
analysis provides an important new approach for the investigation of pathogenic 
mechanisms. 
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The present invention provides greatly improved methods, 
compositions, and apparatus for identifying gene function and for studying 
the regulatory relationship among genes. It is to be understood that the 
above description is intended to be illustrative and not restrictive. Many 
5 variations of the invention will be apparent to those of skill in the art upon 

reviewing the above description. By way of example, the invention has 
been described primarily with reference to the use of a high density 
oligonucleotide array, but it will be readily recognized by those of skill in 
the art that other nucleic acid arrays, other methods of measuring transcript 
10 levels and gene expression monitoring at the protein level could be used. 

The scope of the invention should, therefore, be determined not with 
reference to the above description, but should instead be determined with 
reference to the appended claims, along with the full scope of equivalents to 
which such claims are entitled. 
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We claim: 

1 . A method of determining the stage of disease caused by HCMV infection, 
comprising the step of: 

determining expression levels in a first human cell sample of one or 
5 more genes which are induced or repressed by HCMV, wherein the first 

human cell sample comprises cells of a patient infected with HCMV, 
wherein the first human cell sample consists essentially of HCMV-infected 
cells, wherein the expression levels of the one or more genes correlates with 
stage of disease progression of the HCMV infection. 
10 2. A method of determining the extent of tissue damage caused by HCMV 

infection, comprising the step of: 

determining expression levels in a first human cell sample of one or 
more genes which are induced or repressed by HCMV, wherein the first 
human cell sample comprises cells of a patient infected with HCMV, 
15 wherein the first human cell sample consists essentially of HCMV-infected 

cells, wherein the expression levels of the one or more genes correlates with 
extent of tissue damage caused by the HCMV infection. 
3. A method for screening to identify candidate drugs for preventing disease 
symptoms caused by HCMV, comprising the steps of: 
20 contacting human cells with HCMV and a test agent; 

determining expression levels of one or more genes which are 
induced or repressed by HCMV ; 

identifying a test agent as a candidate drug if the test agent causes 
the human cells to express the one or more genes at a level at which the 
25 human cells express the one or more genes in the absence of HCMV. 



53 



WO 00/011218 



PCT/US99/18772 



4. The method of claim 1, 2, or 3 wherein the one or more genes are selected 
from the group consisting of those genes identified in Table 1 as induced or 
repressed by HCMV. 

5. The method of claim 1, 2, or 3 wherein the one or more genes are induced 
or repressed to a level which is at least two-fold different than the level of 
expression in the absence of HCMV. 

6. The method of claim 1 , 2, or 3 wherein the one or more genes are induced 
or repressed to a level which is at least four-fold different than the level of 
expression in the absence of HCMV. 

7. The method of claim 1 , 2, or 3 wherein the one or more genes are induced 
or repressed to a level which is at least eight-fold different than the level of 
expression in the absence of HCMV. 

8. The method of claim 1 , 2, or 3 wherein the one or more genes are induced 
or repressed to a level which is at least ten-fold different than the level of 
expression in the absence of HCMV. 

9. The method of claim 1 , 2, or 3 wherein the one or more genes are induced 
or repressed to a level which is at least fifteen-fold different than the level of 
expression in the absence of HCMV. 

10. The method of claim 1, 2, or 3 in which the step of detenniiiing is 
performed by measuring amounts of mRNA expressed by the one or more 
genes. 

11. The method of claim 1, 2, or 3 in which the step of detenniiiing is 
performed by measuring amounts of protein expressed by the one or more 
genes.. 

12. The method of claim 1 , 2, or 3 in which the step of determining is 
performed using an array of oligonucleotides. 

13. The method of claim 1, 2, or 3 in which the step of detenmning us 
performed using serial analysis of gene expression. 

14. The method of claim 1 , 2, or 3 in which the step of detenmning is 
performed using hybridization of nucleic acids on a solid support. 
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15. The method of claim 1, 2, or 3 in which the step of determining is 
performed using cDNA which is made using mRNA collected from the 
human cells as a template. 

16. The method of claim 1, 2, or 3 in which a fluorescent label is used to 
determine expression levels. 

1 7. The method of claim 1 , 2, or 3 in which the one or more genes comprise 
HLA-E. 

1 8. The method of claim 1 , 2, or 3 in which the one or more genes comprise 
Ro/SSA. 

19. The method of claim 1, 2, or 3 in which the one or more genes comprise 
lipocortin-L 

20. The method of claim 1, 2, or 3 in which the one or more genes comprise 
cPLA2. 

2 1 . The method of claim 1 , 2, or 3 in which the one or more genes comprise 
COX-2. 

22. The method of claim 1 , 2, or 3 in which the one or more genes 
comprise thrombospondin-L 

23. The method of claim 1, 2, or 3 in which the one or more genes comprise 
MTTF. 

24. The method of claim 1, 2, or 3 in which expression levels of at least 2 
genes are determined. 

25. The method of claim 1, 2, or 3 in which expression levels of at least 5 
genes are determined. 

26. The method of claim 1, 2, or 3 in which expression levels of at least 10 
genes are determined. 

27. The method of claim 1, 2, or 3 in which expression levels of at least 15 
genes are determined. 

28. The method of claim 1, 2, or 3 in which expression levels of at least 20 
genes are determined. 
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29. The method of claim 1 , 2, or 3 in which expression levels of at least 25 
genes are determined. 

30. The method of claim 1, 2, or 3 in which expression levels of at least 30 
genes are detennined. 

3 1 . The method of claim 1 , 2, or 3 in which expression levels of at least 50 
genes are determined. 

32. The method of claim 1, 2, or 3 in which expression levels of at least 100 
genes are determined. 

33. The method of claim 1, 2, or 3 in which expression levels of at least 200 
genes are determined. 

34. The method of claim 1, 2, or 3 in which expression levels of at least 250 
genes are determined. 

35. The method of claim 1, 2, or 3 in which expression levels of at least 500 
genes are determined. 

36. The method of claim 1, 2, or 3 in which expression levels of at least 1000 
genes are determined. 

37. The method of claim 1, 2, or 3 in which the human cells are fibroblasts. 

38. The method of claim 1, 2, or 3 in which the human cells are lymphocytes. 

39. The method of claim 1, 2, or 3 in which the human cells are epithelial 
cells. 

40. The method of claim 1 , 2, or 3 in which the human cells are lung 
epithelial cells. 

41 . The method of claim 1 , 2, or 3 in which the human cells are neuronal 
cells. 

42. The method of claim 1 or 2 further comprising the step of: 

determining expression levels in a second human cell sample of said- 
one or more genes which are induced or repressed by HCMV, wherein the 
second human cell sample comprises cells of said patient, wherein the 
second human cell sample consists essentially of uninfected cells, wherein 
the first and the second human cell sample comprise the same cell type; 
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comparing the expression levels determined in the first and the 
second human cell samples. 

43. The method of claim 1 or 2 wherein the expression levels determined in the 
first human cell sample are compared to expression levels determined for a 

5 reference sample of uninfected human cells of the same cell type. 

44. The method of claim 1 further comprising the step of: 

correlating the expression levels of the one or more genes with the 
stage of disease progression of the HCMV infection. 

45. The method of claim 2 further comprising the step of: 

10 correlating the expression levels of the one or more genes with the 

extent of tissue damage caused by the HCMV infection. 

46. A method for identifying a subset of genes which are improved targets for 
drug development, comprising the steps of: 

(a) comparing expression levels of at least two genes in two cell 
15 samples, wherein the two cell samples are the same but for differences 

caused by a selected environmental, genetic, disease, or developmental 
agent; 

(b) identifying a set of genes whose expression levels differ between 
the two cell samples; 

20 (c) searching a database to identify an unselected environmental 

agent, gene, disease, or developmental stage previously associated with 
expression or altered expression of individual members of the set of genes; 

(d) identifying a common biological feature between the selected 
environmental, genetic, disease, or developmental difference of step (a) and 

25 the unselected environmental agent, gene, disease, or developmental stage 

identified in step (c), wherein identification of a common biological feature 
between two environmental, genetic, disease, or developmental phenomena 
which both affect expression of a common gene identifies the common gene 
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as being a member of a subset of genes which are improved targets for drug 
development. 

47. A computer-readable medium having computer-executable instructions for 
performing steps comprising: 

(a) comparing expression levels of at least two genes in two 
cell samples, wherein the two cell samples are the same but for differences 
caused by a selected environmental, genetic, disease, or developmental 
agent; 

(b) identifying a set of genes whose expression levels differ between 
the two cell samples; 

(c) searching a database to identify an unselected environmental 
agent, gene, disease, or developmental stage previously associated with 
expression or altered expression of individual members of the set of genes; 

(d) identifying a common biological feature between the selected 
environmental, genetic, disease, or developmental difference of step (a) and 
the unselected environmental agent, gene, disease, or developmental stage 
identified in step (c), wherein identification of a common biological feature 
between two environmental, genetic, disease, or developmental phenomena 
which both affect expression of a common gene identifies the common gene 
as being a member of a subset of genes which are improved targets for drug 
development. 
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1 * fx] As all required additional search fees were timely paid by the applicant; this international search report covers all searchable 



claims. 



2 - [ | As all searchable claims could be searched without effort justifying an additional fee, this Authority did not invite payment 
of any additional fee. 

3- | | As only some of the required additional search fees were timely paid by the applicant; this international search report covers 
only those claims for which fees were paid, specifically claims Nos.: 



4 - J^] No required additional search fees wore timely paid by the applicant Consequently, this international search report is 
restricted to the invention first mentioned in the claims; it is covered by claims Nos.: 



Remark on Protest Q The additional search fees were accompanied by the applicant's protest 

[ X| No protest accompanied the payment of additional search fees. 
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International application No. 
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' B. FIELDS SEARCHED 

Electronic data bases consulted (Name of data base and where practicable terms used): 

EAST. DIALOG: MEDLINE, CONF PAPERS, BIOTECH, EURO, JAPIO 

search terms: cytomegalovirus, CMV, HCMV, DNA array, DMA ohip, computer, HLA E, Ro /SSA , Iipocortin-l f 
cPLA2» cystoUe pholpholipase, COX-2, cyclooxygcnase-2, MTTF, microphthalmia associated transcription, fibroblasts, 
lymphocytes, epithelial, lung, neuron?, thrombospondin-1. 

BOX tt. OBSERVATIONS WHERE UNITY OF INVENTION WAS LACKING 
This ISA found multiple inventions as follows: 

This application contains the following inventions or groups of inventions which are not so linked as to form a single 
inventive concept tinder PCT Rule 13.1. In order for all inventions to be searched, the appropriate additional search fees 
must be paid. 

Croup I, claims, drawn to 1-45, drawn to methods of differentiating disease stage or tissue damage caused by human 
cytomegalovirus (HCMV) and to methods of screening for candidate drugs for prevention of symptomatology. 
Group II, claim 46, drawn to methods for identifying a subset of genes which are targets for drug development 
Group m, claim 47, drawn to a computer-readable medium. 

I The inventions listed as Groups Mil do not relate to a single inventive concept under PCT Rule 13.1 because, under 
PCT Rule 13.2, they lack the same or corresponding special technical features for the following reasons: The special 
technical feature of Group I relates to methods of determining disease state or tissue damage from HCMV infection. The 
method of Group II does not possess this special technical feature because it is drawn to a method having different steps 
and a different purpose than the methods of Group L The computer-readable medium of Group III also does not possess 
the special technical feature of Group I because it also relates to a different method having different steps and for 
different purposes than the methods of Group I. 
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